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MOTIVATION 


High  bending  stiffness  and  strength  to  weight  ratio 

Excellent  thermal  and  sound  insulation 

Increased  durability  under  a  thermo-mechanical  loading 
environment 

Tight  thermal  distortion  tolerances 
Lightweight  in  structure 
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BASIC  ASSUMPTIONS  AND  PRELIMINARIES 


1 .  The  face  sheets  fulfill  the  Love-Kirchoff  assumptions  and  are  thin 
compared  with  the  core. 


2.  The  bonding  between  the  face  sheets  and  the  core  is  assumed  to  be 
perfect. 


3.  The  kinematic  boundary  conditions  at  the  interfaces  between  the  core 
and  the  facings  are  satisfied. 


4.  The  core  is  assumed  to  be  a  weak  orthotropic  transversely 
compressible  core  carrying  only  the  transverse  strains  and  the  normal 
strain. 


5.  The  shock  wave  pressure  is  uniformly  distributed  on  the  front  face  of 
the  sandwich  plate. 
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top  face  sheets 


core 


bottom  face  sheets 


Fig  1b.  An  asymmetric  sandwich  plate  under  blast  loading 
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THEORETICAL  DEVELOPMENTS 
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Note: 

the  Greek  indices  have  the  range  1 , 2,  while  the  Latin  indices  have  the  range  1 , 2,  3  and 
unless  otherwise  stated,  Einstein’s  summation  convention  over  the  repeated  indices  is 
assumed.  Also,  denotes  partial  differentiation  with  respect  to  the  coordinates  ,  while 
superscripts  t  and  b  indicate  the  association  with  the  top  and  bottom  facings  respectively. 


Also, 


a  1 /  t  ,  b\  d  1  /  t  b\ 

ui  =  -  (Mi  +  ui  h  =  -  ( ut  -  u i ) 

represent  the  average  and  the  half  difference  of  the  face  sheet  mid-surface  displacements  while,  the  core 
displacements,  warping  functions  of  the  core. 
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Non-Linear  Strain-Displacement  Relationships 

The  strain-displacement  relationships  given  by  the  Lagrangian  Strain-Displacement 
Relationships  used  in  conjunction  with  the  Von-Karman  assumptions  is  given  in 
indicial  notation  as 


1  , 

yu  ~v\,i  +  2^3,1) 

722  =d2,2  +  2^°3,2)2 

t  /  \2 

733  -  d3,3  +  2^3,3) 

__  1  ,  ,1 

723  “2^2,3  +d3,2)+  2D3,2d3,3 

713  _  2  Dl’3  +d37^+ 2D3’lD3’3 

_  1  ,  ,1 

712  -  2  ^Dl’2  +  Vl’1  ^ +  2  D3’lt>3’2 
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Substitution  of  the  displacement  relationships  gives: 
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In  the  above  expressions,  y^,d)  are  referred  to  as  the  average  and  half  difference  of 
tangential  or  membrane  strains  of  the  top  and  bottom  facings;  while,  ^d)  are  referred  to  as 
the  average  and  half  difference  of  the  bending  strains  of  the  top  and  bottom  facings.  The 
expressions  for  the  membrane  and  bending  strains  are  not  provided  here. 


For  the  core,  the  strain-displacement  relationships  take  the  form 


C  -C  ,  c 

7/3  =  7/3  +  ZKi3 

In  these  expressions,  y\ %  and  ^3  are  the  membrane  and  bending  strains,  respectively.  These 
expressions  are  not  provided  here. 
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Constitutive  Equations 

Both  the  top  and  bottom  face  sheets  are  considered  to  be  constructed  from  unidirectional 
fiber  reinforced  anisotropic  laminated  composites,  the  axes  of  orthotropy  not  necessarily 
being  coincident  with  the  geometrical  axes.  The  stress-strain  relationships  for  each  lamina 
of  the  facings  becomes 


Tu  Q\  i  Q\  2  Q\  6  yn 
r22  >=  022  026  '  ^22 
t\2  \  [Sym  066  J  l27i: 


^22 

^12, 


Where,  Qtj  for  i,  j  =  (1 , 2,  6)  are  the  Transformed  plane-stress  reduced  stiffness  measures. 


The  stress-strain  relationships  for  the  orthotropic  core  with  the  geometrical  and  material 
axes  coincident  are  expressed  as 
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Hamilton’s  Variational  Principle 


\,'(SU-SW-ST)dt  =  0 
% 


U  =  strain  energy, 

W  =  represent  the  work  done  by  external  forces 
T  =  represent  the  kinetic  energy 
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Where  Ty  are  the  tensorial  components  of  the  second  Piola-Kirchoff  stress  tensor,  while 
A  is  attributed  to  the  area  of  the  sandwich  plate. 


<5 W  =  $  (?3 (X[,X2 , U&J3  +  ?3  (Xi ,x2 , U&J3  -  2Cti>tiSvti  -  2C%Svi  - 2CbvbSob  }w 

Where  qt(xl,x2,t )  denotes  the  transverse  pressure  loading  from  a  spherical  air-blast 
and  C  is  the  structural  damping  coefficient  per  unit  area  of  the  plate. 


£  dTdt  =  £  L  -  (C/2-^  />%c*3  *3  +f /2+,/  aK  *3 


y 


Where  pc  and  P/ ,  p/  are  the  mass  densities  of  the  core  and  the  top  and  bottom  face 
sheets,  respectively,  v  and  denotes  the  transverse  acceleration. 
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Eq.  (8) 


Where,  the  global  stress  resultants  and  stress  couples  are  defined  as 

+Nb0ji){M>0ji  +Mbaf ,)} 

{NdapMdap)=\ f - N% )  (M^  - M% )} 


UNCLASSIFIED  TECHNOLOGY  DRIVEN.  WARFIGHTER  FOCUSED. 


Physics-Based  Modeling  in  Design  & 
Development  for  U.S.  Defense  Conference 


UNCLASSIFIED 


Where  the  local  stress  resultants  and  stress  couples  are  given  as: 
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0  =  1,  2,  3) 
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Boundary  Conditions 


For  the  case  of  simply  supported  boundary  conditions,  the  boundary  conditions  become: 
Along  the  edges  xn  =  (0,  Ln) 


Kn  =  Kn  =  K  =  K  =  M“n  =  Mi  =  tt?  =4=0 


n  and  t  are  the  normal  and  tangential  directions  to  the  boundary.  When  n  =  1,  t  =  2 
and  when  n  =  2,  t  =  1 
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Solution  Mthodology 
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Special  Case:  Symmetric  orthotropic  single  layer  facings 

In  fulfillment  of  the  geometric  boundary  conditions,  a  suitable  representation  for  u% ,  and  u f 
is  given  by: 


«3  =  M'mnWsin (  lmX[)sin(  nnx2) 
«3  =M'mn(Osin(Amx1)sin(,a„x2) 

Where,  Xm=mKlLl,  nn=n7zlL2 


The  transverse  explosive  loading  is  represented  as 


q,  (x, ,  X2 , 0  =  q„m  (0  S®(  4*1 )  S®(  A„*2  ). 
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which  implies  through  integration  of  both  sides  over  the  plate  area  that 


4 

^?m(0  =  ^-  L  L  qt(xux2it)sm(  4 x, ) sin( nn x2 )dxx dx2 
L^L2  ju  ju 


Letting, 


qt(xl,x2,t)  =  qt(t)  =  (qS0-q0)[l-(t-ta)/tp]exp[-a(t-ta)/tp] 


And  integrating  gives 


(0  =  1H(0 

Hmn\  '  2 
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The  first  two  Equations  of  Motion  can  be  satisfied  by  a  stress  potential  in  conjunction  with  a 
compatibility  equation  not  provided  here.  Equations  of  Motion  (3)  through  (6)  can  be  shown 
to  be  satisfied  by  expressing  these  equations  of  motion  in  terms  of  displacements  and 
assuming  appropriate  functional  forms  in  terms  of  unknown  constant  coefficients  and  the 
amplitudes  as  a  function  of  time.  The  unknown  constants  are  determined  by  substitution 
and  comparing  coefficients. 

At  this  point  the  Extended-Galerkin  Method  is  utilized  by  retaining  the  last  two  Equations  of 
Motion  within  the  energy  functional  and  carring  out  the  indicated  integrations  results  in  two 
nonlinear  coupled  second  order  ordinary  differential  equations  in  terms  of  the  modal 
amplitudes.  These  are  given  as: 


ml  ™mn  +  C™mn  +  ClOwmn  +  C 
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The  coefficients  C10  -C12,C30,C01  -C03,C20,C21  are  expressions  which  depend  on  the 
material  and  geometrical  properties  of  the  structure. 


These  two  governing  differential  equations  are  then  solved  using  the  4th  Order  runge-Kutta 
Method. 
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Blast  Loading 

For  a  free  in-air  spherical  air  burst,  the  pressure  profile  over  time  is  given 
in  figure  2  as 


Fig  2.  Incident  Profile  of  a  blast  wave 
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The  wave  form  shown  in  figure  4  is  given  by  an  expression  known  as 
The  Friedlander  equation  and  is  give  as 


m = (pso  -  p0) 


f  .  .  \  -a 

“I  _  f  ~  *3 
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t-t  ^ 

1  l_a 

lp  j 


P  J 


Where, 


1772 

Z3 


114  108 

— „  H - 


=>  Peak  Overpressure  over  ambient 


Z=  Ft/ l/l/1 3  ^  scaled  distance 


R  is  the  Standoff  Distance 
<  1 N  is  the  equivalent  weight  of  charge 

of  TNT  in  terms  of  kilograms 
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P0  is  the  ambient  pressure 
ta  is  the  time  of  arrival 

tp  is  the  positive  phase  duration  of  the  blast  wave 
f  is  the  time 

For  conditions  of  STP  at  sea  level,  the  time  of  arrival  and  the  positive 
phase  duration  can  be  determined  from 


Arrival  time  or  positive  phase  duration 


'■Mj 


Cube  root  scaling 


Arrival  time  or  positive  phase  duration  for  a  reference  explosion  of  charge  weight, 


l/l/-, 


It  should  be  noted  that  the  standoff  distances  are  themselves  scaled 
According  to  the  cube  root  law 
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Results-  Validation 

To  validate  the  present  approach,  the  dynamic  response  of  a  simply  supported  plate 
impacted  by  a  uniform  pressure  pulse  was  chosen  from  R.S.  Alwar  et  Al. 


Time  (sec) 

Fig  3.  The  nondimensional  global  deflection-time  response  of 
a  simply  supported  sandwich  plate  impacted  by  a  uniform 
pressure  pulse 
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Fig.  4  The  effect  of  the  transverse  modulus  of  the  core  on  the  global 
response  of  a  sandwich  plate  with  orthotropic  facings. 
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Fig.  5  The  counterpart  of  Fig.  4  for  the  wrinkling  response  of  a 
sandwich  plate. 
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Fig.  6The  effect  of  the  rate-of-decay  parameter  on  the  global 
response  of  a  sandwich  plate  with  orthotropic  facings. 
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Fig.  7  The  counterpart  of  Fig.  6  for  the  wrinkling  response. 
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Fig.  8  The  effect  of  the  core  thickness  on  the  global  deflection-time 
history  of  a  sandwich  plate  with  orthotropic  facings. 
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Fig.  9The  counterpart  of  Fig.  8  for  the  wrinkling  response. 
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Fig.  10  The  effect  of  the  stacking  sequence  of  the  facings  on 
the  global  response  of  a  sandwich  plate. 
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Fig.  11  The  counterpart  of  Fig.  10  for  the  wrinkling 
response. 
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Fig.  12  The  effect  of  the  core  shear  modulus  ratio  on  the  deflection-time 
history  of  cross-ply  laminated  sandwich  plate. 
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Fig.  13  The  counterpart  of  Fig.  1 2  for  the  wrinkling  response. 
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Concluding  Remarks 


The  governing  theory  of  asymmetric  sandwich  plates  with  a  first-order  compressible  core 
impacted  by  a  Friedlander-type  of  blast  has  been  presented  and  simplified  for  the  case  of 
symmetric  cross-ply  and  single-layered  orthotropic  facings.  In  all  cases,  it  was  mentioned 
that  all  four  edges  are  simply  supported  and  freely  movable.  Results  were  then  presented 
for  this  simplified  case  and  validated  against  results  found  in  the  literature  from 
R.  S.  Alwar  et  al.  It  was  found  that  for  the  incompressible  core  case  that  there  was  close 
agreement  among  the  results.  In  regards  to  the  compressible  core  case,  no  appropriate 
results  have  been  found  in  the  literature  for  the  theory  presented  in  this  paper  for  the 
simply  supported  case  with  all  edges  freely  movable.  The  effect  of  a  number  of  important 
geometrical  and  material  parameters  were  analyzed  with  conclusions  drawn.  Some  of  the 
important  conclusions  were  that  wrinkling  response  seems  to  be  diminished  as  the  young’s 
modulus  of  the  core  is  increased.  The  same  is  the  case  for  larger  rates  of  decay.  Also,  for 
thicker  cores,  both  the  global  and  wrinkling  responses  are  less  severe.  It  was  also  revealed 
that  the  compressibility  of  the  core  has  only  a  marginal  effect  upon  the  global  response  of 
the  sandwich  plate.  Finally,  the  cross-ply  type  layup  when  compared  with  single-layered 
facings  seemed  to  have  a  large  effect  on  the  global  response  and  less  effect  on  the 
wrinkling  response. 

One  should  keep  in  mind  that  both  the  stress  and  strain  profiles  should  be  determined  to 
determine  possible  failure  of  the  structure. 
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